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We introduce a new class of mesoscopic heat engines consisting of a tunnel junction coupled to a 
linear thermal bath. Work is produced by transporting electrons up against a voltage bias like in 
ordinary thermoelectrics but heat is transferred by microwave photons, allowing the heat bath to 
be widely separated from the electron system. A simple and generic formalism capable of treating 
a variety of different types of junctions and environments is presented. We identify the systems 
and conditions required for maximal efficiency and maximal power. Very high figures of merit and 
efficiencies above the Curzon-Ahlborn limit are achievable. 
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I. INTRODUCTION 

Conversion of heat to work, and, in particular, recov- 
ery of waste heat produced by electronic components is a 
problem of great and ever increasing importance. Solid- 
state thermoelectric systemsi are ideally suited for this 
purpose since they are easily integrated with the rest of 
the circuitry on the microchip. Mesoscopic heat transfer 
devices are a promising class of thermoelectric systems 
due to their ease fabrication, control, and measurement;^ 
and because sharp features in the energy spectrum, a 
requirement for efficient operation,— are readily avail- 
able. Due to their small size they can also be used 
to study foundational issues, such as the importance of 
fluctuations'* and the fundamental limits of heat engine 
performance.^ 

Arguably the simplest mesoscopic heat engine consists 
of a single-level quantum dot placed between two metallic 
leads held at different temperatures and voltages.-. Posi- 
tioning the dot level far enough from the Fermi levels of 
the leads enables the electrons to flow against the voltage 
bias while carrying heat from hot to cold. Thermoelec- 
tric properties of weakly coupled quantum dots have also 
been experimentally studiedj^"— This type of device op- 
erates as a heat engine by generating electrical current 
and transferring heat between the same two reservoirs. A 
recent modificationiiJi of this scheme makes the charge 
and heat currents flow along different pathways by intro- 
ducing a third reservoir: charge is transported between 
two reservoirs at the same temperature, while the third 
reservoir, at a different temperature, is Coulomb-coupled 
to the transport electrons and supplies the thermal fluc- 
tuations driving the heat engine. 

Here we introduce a new type of mesoscopic heat en- 
gine, sketched in Fig. [Ha). It consists of one or more 
quantum dots or metallic islands between two electronic 
leads at the same temperature T but with a voltage bias 
V . The leads are connected to an external circuit with 
temperature Te and impedance Z{lj). If the tunnel cou- 
pling for one of the junctions is small enough, a tunneling 
electron will exchange energy with the electromagnetic 



(a) 




(c) 



A 



FIG. 1. (color online), (a) Schematic of a single-electron heat 
engine coupled to an external environment with impedance 
Z{uj). Electrons tunneling through the left junction exchange 
energy with the heat bath, enabling net current against the 
voltage bias, (b), (c) Energy- level diagrams of junction sys- 
tems with one or two dots between the leads. One of the 
junctions, marked with a photon symbol, is coupled to the 
external bath. The photon-assisted tunneling rates are r±. 



environment, and with a proper choice of parameters it 
is then possible to achieve a net current against the volt- 
age bias. This system is similar in spirit to the devices of 
Refs. [nl and [l3 but instead of a direct Coulomb interac- 
tion between the electrons in the transport system and in 
the heat bath, in our proposal the coupling is mediated 
by microwave photons. Photonic heat conduction in elec- 
tronic circuits has been investigated in recent years both 
theoreticalljiH"— and experimentally ; while previ- 
ously studied single-electron devices for thermal appli- 
cations include a thermometer— i^i^ a cooler—, and a heat 
diodei^ 

A significant advantage of the photonic coupling is that 
the two parts of the circuit with different temperatures 
can be located arbitrarily far apart. One can imagine a 
scenario where the external circuit is a relatively large 
device which performs some useful function but at the 
same time produces excess heat. Our engine can recover 
a part of this heat and feed it back to the main device as 
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electrical power. 

We show that in an optimal configuration the electrical 
current in different types of junction systems is given 
by a single concise formula, Eq. ([8|). Then a junction 
between two quantum dots is shown to be ideal in terms 
of efficiency while maximal power can be provided by 
a junction between a quantum dot and a metal. The 
effect of the environment spectrum is seen to be relatively 
minor. 

The rest of the paper is organized as follows. In Sec.HIl 
we present the theory for electron tunneling in a linear ex- 
ternal environment. In Sec. lIIII we consider bath-assisted 
transport in junction arrays and conclude that in an ap- 
propriate limit the current is given by Eq. for a variety 
of different systems. The heat engine characteristics of 
different junction types are studied in detail in Sec. IIV[ 
with a particular emphasis on performance under max- 
imum power conditions. In Sec. |V] we summarize the 
results and consider experimental prospects. 



II. COUPLING TO THE ENVIRONMENT 

The exchange of energy between the tunneling electron 
and the external environment is treated with the so-called 
P{E) theory. A thorough account of this formalism is 
given in Refs. 23i and ^24,, and in this Section we only 
present the results relevant for the present study. The 
fundamental assumption is that the coupling between 
two electron systems is weak enough so that tunneling 
between them can be adequately described by the Fermi 
golden rule and that the environment relaxation is much 
faster than the tunneling rate. Then within P{E) the- 
ory the tunneling rate Ti^j between the electron systems 
i and j is given by the Fermi golden rule formula with 
the energy-conserving S function replaced by P{E), the 
probability density to exchange energy E with the envi- 
ronment. Thus we have™ (with h = k — e = 1) 
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£j. For the electron density Pi{£) 
we consider two cases, a single-level quantum dot with 
Pi{e) = 5(e), and a metal with pi(e) = Vif{s) where Vi 
is the density of states and /(e) is the Fermi function. 
Similarly the hole density is Pi{£) = 5{e) for dots and 
Pi(e) = Vi[l — /(e)] for metals. Note that in the quan- 
tum dot case we assume that dot i is occupied and dot 
j is empty, otherwise the rate would vanish. The tun- 
neling matrix element between the initial and final states 
is t, which is taken here to be energy independent. The 
Fermi level of a metal or the single level of a dot is p,i 
with possible Coulomb charging energies absorbed into 
it. The heat current emitted by the environment during 
the tunneling process is obtained from the rate formula 



by weighing the integral with Sj — Ei = —£ij'. 
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The P{E) function for an electromagnetic environment 
can be determined by a circuit theory analysis of the 
system. If the junction, which itself has some capaci- 
tance C, is coupled to an environment with impedance 
Z{uj), the total impedance Zt over the junction is C and 
zlJ) in parallel, that is, Zt(w) = [iluC + I / Z {ijj)]-^ . We 
will omit the rather complicated general expression for 
P{E). For our purposes it is sufficient to note that it 
only depends on the environment temperature Te and 
the real part of Zt{uj). Since P{E) is a probability den- 
sity its integral is normalized to unity, and additionally 
the detailed balance for the environment requires that^^ 
P{-E) = e-^/^''P{E). 

We will consider two simple and prototypical environ- 
ments, a harmonic oscillator and an ohmic resistor. The 
form of the P{E) function for an oscillator with frequency 
Wo can be argued rather intuitively. When an electron 
tunnels through the capacitance C, the expectation value 
of the increase in the energy of the electromagnetic field is 
equal to the charging energy Eq = e^/2C, corresponding 
to Ec/ioq emitted photons at Te = 0. At finite temper- 
atures the electron can also absorb energy from the field, 
and the expected numbers of absorbed and emitted pho- 
tons are^ Aa = and Ae = ^{N + 1), respectively, 
with N = (e™o/^^-l)^^ The P{E) function can then be 
written as a double Poisson process for absorbing m > 
photons and emitting m -I- fc > photons, resulting in a 
net emission of E 
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We consider two special cases of this equation, namely 
the limits of large and small wq ■ When wq ^ Te , we have 

N K e and for each k the term with m = max{0, —k} 
dominates the sum while others are exponentially sup- 
pressed. This results in 



PiE) 



\ - e 

fcez 



Ecx\k\ 



\k\\ 



S{E - kujo) (4) 



where 6 is the unit step function. 

In the opposite limit, ujq <C Te, we have ^ ~ 
and the expected photon numbers Aq « EcTe/i^q and 
Ae = Aa + Ec /ojQ become large so that the two Poisson 
processes in Eq. ([3]) can be replaced by normal distribu- 
tions with mean and variance \a/e- The probability to 
emit a net energy i5 is a convolution of these two Gaus- 
sians, itself also a Gaussian;^ 
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where cP' = 2EcTe- An important fact about this equa- 
tion is that it also describes a highly resistive ohmic 
environment, which can be seen by noting that a har- 
monic oscillator with Ziuj) = iujL has 'RGZt{uj) = 
^[d{u! — u!o) + S{uj + u!o)], while a resistor with Z{uj) ~ R 
has ReZt(w) = R/[l + {wRC)^]. Taking the small- 
frequency limit of the oscillator, — 0, and the large- 
impedance limit of the resistor, i? — cx), will both result 
in ReZf(tj) — !• ^S{uj). Therefore both of these systems 
will also have the same P{E) given by Eq. ([SJ, known 
as the high-impedance environment. The spectrum of 
a low-resistance ohmic environment, on the other hand, 
consists of a broadened delta peak at the origin,-22 ap- 
proximating vanishing environmental coupling and will 
therefore not be considered any further. Equations (|4]) 
and (O, representing discrete and continuous environ- 
ments, will be used to investigate the effect of the bath 
spectrum shape on the performance of the heat engine. 

III. CURRENTS IN JUNCTION SYSTEMS 

Since the aim of the present paper is to study ther- 
moelectric power generation in tunnel junctions interact- 
ing with a heat bath through electromagnetic fields, it 
is natural to start with the simplest possible system of 
this type: a single junction between two noninteracting 
metallic leads coupled to an environment with arbitrary 
P{E). The tunneling rate for the positive direction, from 
left (L) to right (i?), is r+ and the rate in the opposite 
direction is r_. With the voltage bias /i_r — /iL = V, 
Eq. ID) gives the current as 

/ = r+ - r_ (6) 

^-i j dsLdeR f{eL)[.f{V - er) - f{~V ~ £i?)]P(e) 

where 7 = 2TT\t\'^i'Li^R and e = El ~ £r- We use a con- 
vention where a positive voltage V means that electrons 
have a higher energy in the right lead and a positive cur- 
rent / means that electrons travel, on the average, from 
left to right. Therefore if / and V have the same sign, 
that is, the power P = IV is positive, heat is converted 
to electrical energy and the device operates as a heat en- 
gine. For the present case the bracketed term in Eq. © 
has a sign that is opposite to the sign of V and therefore 
P < 0. This simple junction cannot produce thermoelec- 
tric power. One way to understand this failure is to think 
of the system as a Brownian motor— which is driven 
by thermal fluctuations due to the external environment. 
However, these fluctuations do not intrinsically have any 
preferred direction and therefore the noise must be recti- 
fied by the electron transport system in order to generate 
net power 1^ Rectification requires a nonlinear current- 
voltage characteristic but the simple junction (without 
the environment coupling) is linear. Therefore a non- 
linearity must be introduced into the transport system, 
and here we do it by adding a quantum dot or a metal- 
lic island between the two leads. We further assume that 



the Coulomb repulsion within the central dot is so strong 
that it can only be empty (with probability Pq) or singly 
occupied (with probability Pi). The Fermi level of the 
metallic island or the single-particle level of the quantum 
dot is at energy E, with Coulomb energy included; see 
Fig. [Ifb). The tunneling rates through the left and right 
junction are T± and F^i, respectively. Now, instead 
of Eq. the current is given by / = PqTj^ — PiF_. 
The probabilities can be solved from the master equa- 
tion Po = Pi(r- + Tr+) - Po(F+ +Tr^) in the steady 
state, Po = 0. This yields the current 

r+ + Fh+ + r_ + Ffi_ ^ ^ 

If the two junctions are identical, symmetry of the system 
implies that there can be no thermally induced current 
for 0, and a calculation with Eq. ([T]) shows that there 
can be no power production for any V . For instance, in 
the case of a quantum dot the numerator of Eq. ([T]) is pro- 
portional to / dede' P{e)P{e')[f{e + E)f[e' - E + V)- 
f{e+E-V)f{e' -E)]. For y > the first and the second 
Fermi function are smaller that the third and the fourth, 
respectively, implying / < 0. To arrive at the preceed- 
ing expression the P{E) functions for the two junctions 
must be identical, and therefore unequal environmental 
couplings are required for power production. 

Since both junctions are connected to the same exter- 
nal circuit, the only way they can have different P{E) 
functions is that the junction capacitances, Cl and Cr, 
arc unequal. In that case circuit theory analysis^"^ shows 
that the total impedance over junction i is KeZt^i^cj) = 
nfRc Zt{uj), where ki = C/Ci and Zt{uj) = [iuC + 
l/Z(a;)]-i, with 1/C = 1/Cl + 1/Cr. In other words, 
the only difference to a single junction is that the ex- 
ternal impedance seen by junction i is effectively re- 
duced by a factor k|. Thus a maximal difference be- 
tween the two P{E) functions can be obtained when the 
capactances have different orders of magnitude. If, for 
example, Cl ^ Cr, we have C « Cl, « 1, and 
kr « Cl/Cr. In this limit the left junction is coupled 
to the environment with energy Ec = e^/2CL while the 
coupling of the right junction is suppressed by a factor 
of (Cl/Cr)"^. Note however that the charging energy is 
e^/2Cfl which must still be large enough to prevent a 
multiple occupation of the island. To simplify Eq. ([7]) 
we note that since the right junction is effectively de- 
coupled from the environment, detailed balance requires 
Tr^ = F7?_e*^'^~^-'/"^. This relation can, of course, be 
also derived from Eq. ([T]) with P{E) = S{E). Further- 
more, a junction with a small capacitance is usually also 
weakly transmitting and therefore it is consistent to as- 
sume that Cl ^ Cr implies T± <C Tr±. Then we obtain 

I ^ fiV -E)T+- fiE~V)T^ (8) 

This is the main equation of the present paper. Below 
we argue that it applies also to more generic junction 
systems, and its heat engine characteristics are examined 
in detail in the next section. 
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Let us then turn to the case of two metal islands or 
quantum dots in series between the leads, as depicted in 
Fig. (He). We again assume a strong enough Coulomb 
repulsion so that the double dot can be either empty 
(probability Pq), or there can be one electron in the left 
dot (probability Pl) or in the right dot (probability Pr). 
The tunneling rates through the left, center, and right 
junction are Tl±, T±, and Tfi±, respectively, and the 
current is given by / = Pl^+ — PrT^. Master equation 
for the occupations is 

Po - -(ri+ + rR^)Po + Tl-Pl + Tb+Pr, 

Pl = Tl+Po - {Tl- + r+)PL + T^Pr (9) 
Pr = Tfl.Po + T+Pl - {rR+ + r_)Ffl 
and the steady-state solution gives the current as 

I ^{TL+r+TR+-TL-r^TR^)/f' (10) 
with f 2 rL+rR+ + Tl-Tr^ + rL-rR+ + r+(ri+ + 

Tr+ +rR-)+ T-{Tl+ + Tl- +Tr-). Similarly to the 
two-junction case, we now assume that the central junc- 
tion has a small capacitance and a small transmittance 
compared to the other junctions. Thus r± ^ Tl±.r±, 
and since the left and right junctions are effectively de- 
coupled from the environment, detailed balance implies 
TL-/rL+ = e^-/^ and Tr+/Tr. = e^^-"^)/^. The 
expression for the current is then simplified to 

1 + e(-E-V)/T _^ ^{En-V)/T ^'^'^) 

The last term of the denominator is vanishingly small 
compared to the first two if —{Er ~~V) ^ T and —El ^ 
T. In this limit the expression for the current is again 
reduced to Eq. ([5]). Physically this limit means that the 
dot levels are much below the Fermi levels of the leads 
and thus either of the two dots is always occupied {Pq — 
0), leading effectively to a two-state system similar to 
the one-dot case. Clearly this argument could also be 
extended to a larger number of dots if necessary. 

We have now seen that Eq. ([8]) is the fundamental ex- 
pression for thermoelectric current in the optimal limit 
when only one junction is exchanging heat with the ex- 
ternal bath and when all "idle" time spent on processes 
with no environment coupling can be neglected. The gen- 
erality of Eq. ([8]) can also be intuitively seen as follows. 
The system is always in one of two possible states: there 
is an electron ready to tunnel either from left to right 
or from right to left through the bath-coupled junction. 
The probabilities of these two states are P_|_ and P^, 
and the average current through the junction is there- 
fore / = P+T^ — P_r_. On the left side of the junc- 
tion there is a metal with Fermi level El^ or a quan- 
tum dot with a level at El strongly coupled to a metal. 
In either case, the probabilty that the left side has an 
electron ready for tunneling is proportional to /{El)- 
Similarly the right side has a level at Er — V, and it 
is ready to receive the tunneling electron with a prob- 
ability proportional to 1 — /{Er — V) = f{V — Er). 



Thus we have P+ oc f{EL)f{V — Er) and analogously 
P- cx f{-EL)f{ER-V). Normalizing with P+ + P^ = 1 
gives Pq = f{V — E) and Pi = f{E — V), and we arrive 
at Eq. The existence of these two states is due to the 
electron-electron interaction. The noninteracting case of 
Eq. (IH) is obtained when P+ = P_ = 1, that is, when 
there is only one single state allowing electrons to tunnel 
at any time in either direction. 

The heat current J emitted by the environment can be 
calculated by following exactly the same steps as for the 
electrical current, and the result corresponding to Eq. (|8]) 
is 

J = f(V- E)J+ + f{E - V)J- (12) 

where J+ and J„ are the energy currents absorbed by 
the electron tunneling to the right and to the left, re- 
spectively, as obtained from Eq. ([2]). The sign difference 
between Eqs. ([8|) and (fT2|) is due to the fact that J± 
already contain the direction of the heat flow. 

IV. HEAT ENGINE CHARACTERISTICS 

As we have argued, in the optimal limit an array of 
one or more quantum dots or metal islands can be seen 
as a two-state system obeying Eqs. ^ and ([T2|) . and thus 
when evaluating the thermoelectric power generation in 
these devices it is only necessary to consider the single 
junction that is coupled to the electromagnetic environ- 
ment. Some generic remarks about Eq. ([8|) can be made 
without any explicit model for the junction. First note 
that a mirror reflection of a solution with {E, V, I) pro- 
duces another solution with {—E, —V, —I) and therefore 
it is sufficient to consider the case E > 0. Next, the I{V) 
curve of Eq. ([5]) has a simple overall structure. Since the 
tunneling rates T± depend only on the level difference 
E over the junction but not on the voltage bias V, the 
current depends on V only through the Fermi functions. 
Thus we see that / approaches asymptotically F-i- and 
— r_ for large negative and positive values of V, respec- 
tively. Since I{V) decreases monotonically, at some point 
V = Vq the current vanishes. For our purposes the most 
important fact is that for voltages between and Vb, / 
and V have the same sign and the device operates as a 
heat engine. From Eq. ([8]) we can solve 

Vo^E-T\og^ (13) 

All the systems studied here have r_ > r+ and therefore 
E is an upper limit for Vq while no lower limit exists. We 
will also see that a hot environment {Te > T) implies 
Vb > and therefore the heat engine operates with pos- 
itive V and / while for a cold environment [Te < T) V 
and / are negative. 

An important characteristic of a heat engine is the ef- 
ficiency r/ — IV/ Jh^ where Jh is the heat current from 
the hot bath. When Te > T we have simply Jh — J 
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but when Te < T the hot bath is the transport system 
and then Jjj — IV — J, that is, the heat taken from the 
electrons is the sum of the produced power IV and the 
heat — J expelled to the cold environment. Thus we have 
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(14) 



The fundamental upper limit for r/ is the Carnot effi- 
ciency r]c — I ~ Tc/Th, where Tc = minjT, Tg;} and 
Th = max{T, Tb}. In this limit transport proceeds re- 
versibly and thus power production is vanishingly small. 
In practice one often wishes to operate the engine with 
parameters that maximize the power instead of efficiency. 
As shown by, among others, Curzon and Ahlboru;^ us- 
ing a particular model for a heat engine it is possible to 
derive a simple upper bound 7]ca = 1 — \/Tc/Th for the 
efficiency at maximum power. However, later it has been 
realized that rjcA is not a universal limit;^ and in fact, 
as we will see below, our engine is also able to exceed it. 

Another efficiency measure, widely used in the con- 
text of thermoelectric power generation, is the figure of 
merit ZT. Even though the separation of heat and charge 
pathways makes the devices covered by our theory quite 
different from typical thermoelectrics, the usual defini- 
tion of ZT can still be straightforwardly applied to the 
present case. In the linear response limit, when V and 
AT = Te — T are much smaller than T, we have^ 



ZT 
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(15) 



where the Seebeck coefficient is S* = 9Vb/9(Ar), elec- 
trical conductance is ct = dl/dV at AT ~ 0, and ther- 
mal conductance is k — dJ/d{AT) at / = 0. Gener- 
ally K should include all forms of heat transfer between 
the reservoirs, but since we are not modeling any para- 
sitic flows they are not included in ZT. A central goal 
in thermoelectrics research is to produce systems with 
ZT > IJ. 

When an electron is transported between the leads 
of a thermoelectric system it always performs the same 
amount of work ±V where the sign depends on the di- 
rection of tunneling. However, the amount of heat trans- 
ferred between the thermal baths can vary from one tun- 
neling event to another, and the spectrum of this en- 
ergy exchange is determined both by the structure of the 
junction and the environment. We consider two limit- 
ing cases for the junction structure. The first is a fully 
energy-selective junction, where each tunneling event is 
accompanied by an exchange of a fixed amount of heat 
±E. A physical realization for this is a junction between 
two quantum dots with sharply defined energy levels sep- 
arated by E. The other extreme is a totally unfiltered 
junction where essentially any amount of heat can be ex- 
changed. The two physical examples that we consider 
are a junction between a quantum dot and a metal, and 
a junction between two metals. 



A. Energy-selective junctions 

A junction between two quantum dots can only ex- 
change a fixed energy E with the external bath during 
the tunneling events. Using Eq. ([T]) together with de- 
tailed balance for the environment yields 
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Similarly the heat flows from Eq. 1^ are J_|_ = E^-^- 
and J- = —ET-, which imply the simple relation J = 
EI. Therefore at all temperatures the ratio of produced 
power and transferred heat is a constant V/E, a situa- 
tion known as strong coupling between particle and heat 
flows. Carnot efficiency can only be achieved by such 
strongly-coupled systems ^"2^12^ and in the present case 
this can be confirmed by noting that the stopping voltage 
from Eq. ^ is 



Vo = E{1 - T/Te) 



(17) 



and substituting V — Vq in Eq. (|T4)) yields rj = rjc- Thus 
Carnot eflticiency implies vanishing current and power. 

In the linear response regime, that is, to first order in 
V and AT, the current from Eqs. ^ and (fTO)) is 



I^ll^L^(EAT-TV) 
1 ^ 



(18) 



The value of ZT can be inferred directly by noting that 
since k in Eq. ([TSl) is evaluated at / = 0, the strong- 
coupling condition implies J — k — and therefore 
ZT = oo, independent of any parameters. This is an- 
other indication that energy-selective junctions can be 
used to construct maximally efficient thermoelectrics. 

In order to investigate the physics beyond linear re- 
sponse and to find out the conditions that maximize the 
power instead of efficiency, we perform a numerical cal- 
culation and for that purpose an explicit expression for 
the environment spectrum is needed. Since an energy- 
selective junction interacts with the environment only 
at a single energy E, the full form of the P{E) func- 
tion does not generally have much significance. There 
is one caveat, however: Fermi golden rule, which P{E) 
theory is based on, is not able to treat transitions be- 
tween discrete states, and therefore our approach fails if 
the environment and both sides of the junction are dis- 
crete. Thus for a junction between two quantum dots we 
only consider the continuum environment, Eq. ([S]). Note, 
however, that even in this case the limit ct — ?> gives a 
delta function which is beyond P{E) theory. 

For each pair of temperatures T and Te, the generated 
power P = IV is numerically maximized with respect to 
the bias V, the level difference E, and the coupling en- 
ergy Ec, and the results are presented in Fig. [5Ja) as a 
function of AT/Tq, where Tq = ^{T + Te) is the aver- 
age temperature. Beyond linear response the behavior of 
the system becomes non-symmetric in AT. The limit- 
ing cases of large temperature differences can, however. 
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FIG. 2. (color online). Maximum power production for (a) a 
junction between two quantum dots coupled to a continuum 
environment, and (b) a junction between a quantum dot and a 
metal coupled to a discrete environment with ljq ~ STq. Top 
panels: The maximum power for given temperatures. Note 
the different units for (a) and (b). Middle panels: efficiency rj 
(full line) and the efficiency compared to the Curzon-Ahlborn 
value, rj/rjcA (dashed), at maximum power. Lower panels: 
values of V (full line), E (dashed), and Ec (dash-dotted) 
that give the maximum power. The temperature difference is 
AT = Te-T and the average temperature is To = ^{T + Te). 



be understood from simple analytical arguments. When 
T (AT ^- 2To), we have F > and therefore the 
current must also be positive. In this limit the Fermi 
functions in Eq. ([8|) become step functions and we must 
have E' — > so that the forward tunneling process re- 
mains while backward tunneling vanishes. On the other 
hand, to maximize power V should be as large as possi- 
ble, and therefore we must have V E. This is also the 
condition for maximum efficiency, since in this case all 
the heat E taken from the hot environment is used for 
work V, implying 77 — 1. This does not violate the sec- 
ond law since a vanishingly small energy E—V emitted to 
a cold bath with T — ^ is enough for the entropy change 
AS* = -~E/Te + {E-V)/T to be positive. Thus when the 
cold bath is approaching zero temperature it is possible 
to have efficiencies close to the Carnot value with finite 
power. Finally, the actual value of the maximum power 
and the values of Ec and E = V are determined by the 
form of P{E). The calculation shows that the optimal 
case for Eq. ([5]) has Ec and therefore tr — >■ 0. As 
discussed above, the environment spectrum now becomes 
a delta peak which is beyond the present theory, leading 
to an unphysical divergence in the power. This prob- 
lem only occurs when the two temperatures have very 
different magnitudes and therefore it is of no practical 
consequence. 



In the opposite limit, Te — )• (AT — > — 2To), we have 
V < Q and and the current must also be negative. The 
requirement that E > Q implies r+/r_ = e^^/^^ -> 
so that tunneling in the positive direction is suppressed. 
This is because the energy-selectivity of the junction 
forces the electron to absorb an energy E from the envi- 
ronment during this process, but in the limit Te — this 
energy is not available. From the second Fermi function 
of Eq. ([H) we then see that the magnitude of the current, 
and thus the power, is maximized when E is as small as 
possible, and therefore the condition for maximum power 
is E 0. This also maximizes the efficiency since only 
a vanishingly small energy E needs to be expelled into 
the cold environment so that in order to do work \V\, 
only the same amount of energy needs to be taken from 
the hot electron system, again resulting in 77 — > 1. The 
remaining parameter values are determined by numerical 
calculations but just as in the previous case, also here the 
maximization with Eq. ([5]) results in cr and in the 
unphysical divergence of power. 

One more indication of the highly efhcient nature of 
the energy-selective junction is the fact that the value of 
77 at maximum power stays rather close to the Curzon- 
Ahlborn limit. The efficiency at maximum power has 
been actively studied in recent years and it has been 
shown^^'^^ that for a system with strongly coupled 
heat and particle flows the efficiency will approach the 
Curzon-Ahlborn value in the limit AT — > 0. Such gen- 
eral results are not known for finite temperature differ- 
ences, and as can be seen from Fig.[5fa), the present sys- 
tem actually has rj > rjcA for a cold environment while 
"qcA is an upper bound in the presence of a hot environ- 
ment. 

In addition to the different measures of efficiency, a 
very important characteristic is the absolute magnitude 
of power achievable with the device. The unit of power 
in Fig. [21a) is 27r|ip, and to find an upper limit for this 
quantity we note that the junction between two quantum 
dots can be seen as a two-state system with level differ- 
ence E and tunnel coupling t. The left and right dots 
must be approximate energy eigenstates, requiring that 
\t\ ^ E. Since E Tq for maximum power operation, we 
end up with \t\ ^ Tq- Thus the output power is strictly 
limited by the operating temperature. 

Finally we remark that this energy-selective system is 
very similar to the one studied in Ref. [ll|. The physical 
implementations are rather different, with the system of 
Ref. |T2 requiring a total of four quantum dots while for 
our device two dots are sufficient, but the structures of 
the energy transfer processes are fundamentally identical, 
leading to very similar performance figures. 



B. Unflltered junctions 

Junctions which do not place any restrictions on the 
amount of energy exchanged between the tunneling elec- 
tron and the environment will be called unfiltered. We 
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mostly concentrate on a system where the other side of 
the junction is a quantum dot and the other side is a 
metal, in which case Eq. ([T]) gives 



r± =ry" deP(e)f{e±E) 



(19) 



where T = 2'r:\t\'^v. This form shows immediately that 
r_ > r+ and thus Vb < E. Linearizing Eqs. dH) and 
(fT2)) then yields 



I = -^-^ {El\T - cTV) 



J 



2^2 



rp2 



{ J AT - ETV) 



(20) 
(21) 




while the corresponding figure of merit from Eq. ()15p is 

-1 

(22) 

where c = F[0], E ^ F[l], J = i^[2], and F\n\ = 
deP(e)e"[/(£:-e)e-T + (-i)"/(£;+e)]. Thesedefini- 
tions show that if PiE) consists of a pair of delta peaks 
at energies ±£^0: then in the limit E = Eq ^ T we 
have ZT oo and therefore this kind of environment 
is able to mimic the effect of an energy-selective junc- 
tion. However, in this limit E, and thus the power IV, 
scales as exp{—E/T). A more detailed analysis shows 
that beyond ZT ~ 1 a linear increase in ZT corre- 
sponds to an exponential suppression of generated power. 
This should be contrasted to the energy-selective junc- 
tion where the infinite ZT is due to the idealized as- 
sumption of perfectly sharp energy levels. In reality the 
levels are broadened and the figure of merit is finite. 
However, when the level width is decreased, power is 
reduced roughly inversely, and not exponentially, with 
30 ^]^yg g^jj energy-selective junction is the only re- 
alistic way of achieving very high efficiencies, at least in 
the linear regime. The fundamental difference between 
the energy filtering provided by the junction and by the 
environment is the fact the former is directional while the 
latter is not: when an electron tunnels to the right in an 
energy-selective junction it must always absorb a photon 
and when tunneling in the other direction it must emit a 
photon, but for an unfiltered junction both emission and 
absorption are possible for either direction. 

After concluding that large efficiencies are not avail- 
able for linear response, we turn to the nonlinear regime 
and conditions for maximum power. Figure [Hb) shows 
the results for a discrete environment as represented by 
Eq. Q with Wo = 5Tq. The figure of merit at maximum 
power is ZT « 0.5. In the limit of large temperature 
differences we see that the case of a hot environment, 
T — >■ or AT — >■ 2To, is very similar to the energy- 
selective junction. Indeed the physics is essentially the 
same: when V E = Eq = wg, only the forward process 
of Eq. ([5]) remains, always absorbing a quantum ujq from 



the hot environment and performing the same amount of 
work, leading to 77 — >■ 1. In the opposite case of a cold 
environment, Te ^ 0, performance of the heat engine is 
dramatically degraded. Since now y < 0, electrons tun- 
neling to the left perform useful work, and for the energy- 
selective system we were able to have high efficiency and 
power because r+ was vanishing. But a detailed inspec- 
tion of Eqs. ^ and reveals that now r_|_ remains 
large because electrons are able to tunnel to the right by 
using the thermal energy of the electron system, without 
absorbing anything from the external bath. The trans- 
port processes which are decoupled from the environment 
produce a considerable leakage current down the voltage 
bias, thus making heat engine performance very poor. 

For practical purposes an ohmic environment is an im- 
portant special case. This is modeled with the high- 
impedance P{E) function from Eq. ([5]), which, as we 
have noted, also represents the single-mode environment 
in the limit 0- Numerical results for power maxi- 

mization are shown in Fig.[3][a). In the linear regime the 
figure of merit is ZT « 0.5. For Te > T the maximum 
power curve is very similar to the case of discrete environ- 
ment, and indeed a calculation shows that for ujq < 5Tq 
the maximum power varies only weakly with loq while for 
higher resonance energies it is exponentially suppressed 
in loq/Te- The efficiency is smaller for the continuum 
environment since the power is generated by elementary 
tunneling processes each performing the same work V 
but absorbing a different amount of heat from the hot 
bath. These elementary processes thus have different ef- 
ficiencies, each limited by the second law, and the over- 
all average efhciency must therefore be lower than the 
Carnot value. 

The case of a cold environment, Te < T, is somewhat 
surprising since both power and efficiency are higher than 
for a discrete environment. The reason for the very poor 
performance of the discrete case was seen to be the tun- 
neling processes without environmental coupling. But for 
the high-impedance P{E), with Te <C Ec, the environ- 
ment spectrum is sharply peaked about Eq, away from 
the origin, and the decoupled processes are suppressed. 
Consequently the maximum power in the Te — >■ limit 
is increased by an order of magnitude. 

As seen from Figs. [Ifb) and[3l^a), for a junction be- 
tween a metal and a quantum dot the magnitude of the 
maximum power depends on the product of F and Tq. 
Since we are assuming that only a single level of the dot 
is accessible, both of these quantities must be small com- 
pared to the dot level spacing but F is not directly con- 
strained by Tq. Therefore in a system with a level spacing 
much larger than the temperature it is possible to have 
TTo » Tl 

Finally we consider another model system for the un- 
filtered junction, namely tunneling between two metals. 
In this case Eq. ([1]) yields 



(23) 



r± = 7 / deP{e)en{e ± E) 
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FIG. 3. (color online). Maximum power production for (a) 
a junction between a quantum dot and a metal coupled to 
a continuum environment, and (b) a junction between two 
metals coupled to a continuum environment. See Fig. [2] for 
an explanation of the different plots. 



where 7 = 27r|^p^'iZ/fl and n{£) is the Bose function at 
temperature T. Comparing this to Eq. we see that 
the quantitative results for the junction between a quan- 
tum dot and a metal can be transferred to the present 
case by replacing F — >• 7 and /(e) — >■ £n(e), while all qual- 
itative arguments remain unchanged. We only present 
the maximum power graph for an ohmic environment in 
Fig.[3jb) and note the considerable similarity to Fig.[3l[a). 
The efficiency is somewhat lower for the present case, 
with ZT sa 0.2. The unit of power in Fig. ^h) is -^T^, 
and the existence of Coulomb blockade requires^ 7 <^ 1. 
Therefore just like in the case of the junction between two 
quantum dots, the operating temperature strictly limits 
the attainable power. 



V. DISCUSSION 

In the previous section we have shown that highly effi- 
cient energy conversion is in practice only available for an 
energy-selective junction. This type of device also per- 
forms equally well for hot and cold environments. Unfil- 
tered junctions should preferably be operated with a hot 
environment and cold electron system, and in this case a 
discrete environment can offer a slightly higher efficiency. 
If, however, this system is used with a cold external bath, 
a continuum environment offers better performance. 

To intuitively understand the behavior of single- 
electron heat engines at large temperature differences, 
one can consider the system in Fig. [Ifb) with the left 
junction coupled to the environment. When the electron 



system is cold, that is, T is small compared to the other 
energies, then the Fermi functions are sharp and electrons 
cannot tunnel up in energy from the right lead to the cen- 
ter, and therefore only the positive direction transport 
processes remain. First an electron tunnels from the left 
lead to the center by absorbing a photon, and then due to 
the strong tunnel-coupling of the second junction it dis- 
charges to the right lead. Thus each tunneling electron 
transfers heat from hot to cold and performs useful work, 
leading to optimal thermoelectric performance. On the 
other hand, for a hot electron system the Fermi func- 
tions are smeared and in general electrons can tunnel in 
both directions without absorbing or emitting photons. 
Such environment-decoupled processes transport electri- 
cal current down the voltage bias, that is, they produce 
Joule heating from work, hence severely degrading the 
engine power and efficiency. Only a fully energy-selective 
junction is able to force the tunneling electrons to al- 
ways interact with the environment, resulting in highly 
efficient energy conversion. 

Typically achieving high efficiency is a major goal in 
heat engine design but for the purposes of waste heat 
recovery the thermal input energy of the device can be 
considered free and abundant, making efficiency an ir- 
relevant quantity. In this case one should instead con- 
centrate on maximizing the power output, as we have 
done in Figs. [H and |31 We showed that for a junction 
between two quantum dots or between two metals the 
output power is limited by the temperature in such a way 
that the units of power in Figs. [IJa) and Ellb) must be 
much smaller than Tq . The same kind of limitation does 
not apply for a junction between a metal and quantum 
dot and thus for fixed temperatures this system has the 
potential for largest power generation. The reason is that 
for a dot-metal junction the coupling strength is limited 
by the dot level spacing which can be much larger than 
To, while the coupling in a dot-dot junction is limited 
by the level difference E which cannot be much larger 
than the temperature or current becomes suppressed. In 
a metal-metal junction the coupling must be much below 
the conductance quantum or otherwise Coulomb block- 
ade is wiped out, leading to a linear junction with no 
thermoelectric effect. If F ~ Tq = 1 -fC, the maximum 
power for all three junction types falls in the femtowatt 
range. This is a typical figure for low-temperature single- 
electron devices 

In principle it is possible to increase the power pro- 
duction by having several engine units in parallel with 
the external impedance. The coupling energy Ec scales 
inversely with the number of parallel devices and there- 
fore the capacitances of the individual junctions should 
be decreased. However, in order to see an actual increase 
in current and power, this change in the capacitances 
should not considerably lower the tunneling rates of the 
junctions. 

Even if one is not concerned with efficiency, heat leaks 
between the thermal baths should be minimized in or- 
der to have a maximally large temperature difference. 
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Phononic heat conduction is a problem with ah ther- 
moelectric systems. Since phonons are not part of our 
model we only note that at very low temperatures the 
electron-phonon coupling becomes weak and this con- 
duction channel can be ignored. On the other hand, two 
heat leak mechanisms particular to the type of devices 
considered here are heat conduction by electrons mov- 
ing between the junction and the external circuit, and 
photonic heat transfer between the external impedance 
and the metallic reservoirs of the junction system. The 
first leakage channel can be eliminated by using super- 
conducting wires, and the second one is suppressed if the 
resistance of the electron reservoirs is much smaller than 
the tunnel resistance of the junctions. 

One should also note that since the generated current 
flows through the external impedance Z(lu), the power P 
and voltage bias V must be related by Z{lu ^ 0) = V"^ /P 
if there are no other voltage sources. For example, from 
Fig. Elb) we then have Z(0) ~ 10^7-1 with 7 < 1 so 
that the external resistance must be several orders of 
magnitude larger than the resistance quantum, which is 
consistent with the assumption of a high-impedance envi- 
ronment. Of course a non-thermoelectric voltage source 



in series with the heat engine can be used to drive the 
current through an even larger load. 

In conclusion, we have investigated the heat engine 
performance of several types of single-electron junctions 
coupled to linear electromagnetic environments. Highest 
thermoelectric performance is obtained when only a sin- 
gle junction is coupled to the external heat bath, and in 
this case a simple formula is able to describe the essen- 
tial dynamics. It was confirmed that an energy-selective 
junction between two quantum dots is capable of highly 
efficient energy conversion, as evidenced by the high ZT 
number and the efficiency at maximum power being close 
to the Curzon-Ahlborn value. However, a junction be- 
tween two metals can reach the same magnitude of max- 
imum power, while highest power generation is available 
with a junction between a metal and a quantum dot. 
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